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Abstract. We study the small scale anisotropy signal expected at the Pierre 
Auger Observatory in the next 1, 5, 10, and 15 years of operation, from sources 
of ultra-high energy (UHE) protons. We numerically propagate UHE protons 
over cosmological distances using an injection spectrum and normalization that 
fits current data up to ~ 10 20 eV. We characterize possible sources of ultra- 
high energy cosmic rays (UHECRs) by their mean density in the local Universe, 
p = 10 _r Mpc -3 , with r between 3 and 6. 

These densities span a wide range of extragalactic sites for UHECR sources, 
from common to rare galaxies or even clusters of galaxies. We simulate 100 
realizations for each model and calculate the two point correlation function for 
events with energies above 4 X 10 19 eV and above 10 20 eV, as specialized to the 
case of the Auger telescope. We find that for r J> 4, Auger should be able 
to detect small scale anisotropics in the near future. Distinguishing between 
different source densities based on cosmic ray data alone will be more challenging 
than detecting a departure from isotropy and is likely to require larger statistics of 
events. Combining the angular distribution studies with the spectral shape around 
the GZK feature will also help distinguish between different source scenarios. 
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1. Introduction 

There are two main open questions that together constitute the long standing mystery 
of the origin of ultra- high energy cosmic rays (UHECRs): 1) What is their acceleration 
site? 2) Does their diffuse spectrum at the Earth possess the so-called GZK feature 
, evidence of their inelastic interaction with the photons of the cosmic microwave 
background (CMB)? In the next years the Pierre Auger Observatory [2] will provide 
us with important hints to the first question and will most likely answer the second 
question. Moreover, it can be expected that this experiment will open the way to 
realistic measurements of the chemical composition of cosmic rays at the highest 
energies, another important piece of the puzzle. 

As the Southern site of the Auger Observatory nears completion, an increase by 
one order of magnitude of the worldwide exposure to the highest energy particles 
is within reach. The energy range between 10 19 eV and 10 20 eV is of particular 
importance as it provides a window for charged particle astronomy. Cosmic rays 
in this energy range suffer smaller deflections by cosmic magnetic fields than their 
lower energy counterparts, although at the lowest end of this range the effects of the 
galactic magnetic field are expected to be important. However one can imagine to be 
potentially able to correct for this effect or even use it to infer the topology of the 
magnetic field of the Galaxy. 

A key element in determining the origin of UHECRs is the detection of large [3] 
and small scale anisotropies in their arrival directions. If, as expected and as suggested 
by some numerical simulations (see, e.g., py), the deflection of cosmic rays of energies 
larger than ~ 4 x 10 19 eV should be small, then to some extent they should point 
back to their source. As large numbers of cosmic rays at such energies are observed, 
the pointing becomes more easily feasible. The net result is the appearance of small 
angle clustering in the arrival directions, and the strength of the effect is related to 
the mean density of sources in the universe. 

In addition, interactions of cosmic rays with the cosmic background radiation 
limit the distance from which cosmic rays with energies above ~ 10 20 eV can reach the 
Earth to ~ 50 Mpc. The matter distribution within a 50 Mpc radius from the Earth 
is rather anisotropic, thus large scale anisotropies in the UHECR arrival direction 
distribution should also be observable. 

An important piece of information, necessary to infer the origin of UHECRs 
is their dominant chemical composition as a function of energy. The present 
observational information about the composition at the highest energies is rather 
poor and does not allow us to draw any firm conclusions. While our calculations are 
carried out for a purely proton composition at the highest energies, an appreciable 
contamination of heavier nuclei may change our conclusions on the implications of 
small angle clustering to some extent. From the theoretical point of view, a proton- 
dominated or a mixed composition also have quite different implications on the 
transition of cosmic rays from a galactic to an extragalactic origin: while a mixed 
composition would imply that such transition takes place at the ankle |SJ [S] , a pure 
proton composition (or very light pollution from heavier nuclei) would imply that the 
transition occurs at lower energies {dip scenario) 

Our conclusions regarding small scale anisotropies should apply more generally 
with the caveat that we assume the best fit injection spectrum for a pure proton 
source, which is softer than the mixed composition case (however, see |S] for the 
effect of magnetic fields in the source region and |§] for the effect of maximum energies 
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depending on the source). Given that the relevant energy range is limited (between 
4 x 10 19 eV and 10 20 eV), small differences in the injection should not have a significant 
effect. 

The unknown magnitude and structure of the magnetic fields outside the Milky 
Way limit our ability to predict the precise energy threshold for charged particle 
astronomy. The role of extragalactic magnetic fields in changing the spectrum and 
clustering of UHECRs is also badly constrained as different simulations give different 
estimates for the magnitude and spatial structure of these fields (see, e.g., The 
combined effect of a mixed composition and significant extragalactic fields may move 
the threshold for detecting anisotropies to larger energies (see, e.g., QD). Here, we 
assume that extragalactic magnetic fields can be neglected for particles with energies 
above 4x 10 19 eV. This assumption holds if magnetic fields in the extragalactic medium 
are less than ~ 0.1 nG with a reversal scale of ~ 1 Mpc and the small scale anisotropies 
are evaluated on angles of ~ 1° \V2i I13|. This magnitude field is compatible with 
observational bounds and detailed numerical simulations such as 0] (however, see 
[TU] for different numerical estimates). UHECR observations will eventually allow us 
to measure the strength of the extragalactic magnetic field and get information about 
its structure (see, e.g., ^JE])- In the mean time, searches for small scale anisotropies 
should focus on energy thresholds from 4 x 10 19 eV to 10 20 eV. 

The mean density of UHECR sources can be determined once small scale 
anisotropies (SSAs) in the arrival directions are observed. Thus far the only 
experiment to report departures from isotropy is AG AS A 17. The statistical 
significance of the clustering has been challenged ^Hj as it depends on the angular 
scale chosen for binning the data. Assuming the AGASA data, the number density 
of sources was estimated to range from ~ 1CT 6 Mpc" 3 to ~ 1CT 4 Mpc" 3 with 
large uncertainties EH 1191 I20| . In Refs. |121 I13| . full numerical simulations 
of the propagation were performed which account for the statistical errors in the 
energy determination and the AGASA acceptance. The first study concluded that 
the AGASA small scale anisotropies indicated a density of sources ~ 10 _5 Mpc -3 |12| . 
However, when taken together with the observed AGASA spectrum, the SSAs and 
the GZK feature become inconsistent 13 . 

Here we study the small scale anisotropy signal expected at the Pierre Auger 
Observatory in the next 1, 5, 10, and 15 years. We numerically propagate ultra high 
energy protons over cosmological distances and characterize possible sources by their 
mean density in the local Universe, p — 10~ r Mpc -3 , with r between 3 and 6. These 
source densities span the relevant range of extragalactic sites for UHECR sources, from 
galaxies to clusters of galaxies. In §2 we review our numerical approach and present 
the SSA results of 100 realizations of each model. We show the two point correlation 
functions for events of energies above 4xl0 19 eV and above 10 20 eV, list the expected 
number of doublets, and discuss the expected spectrum for different models. In §3, 
we discuss the implications of our results and conclude. 

2. Small Scale Anisotropies and Spectra at Auger 

The propagation of UHECRs is simulated using the Monte-Carlo code described in 
\'21\ . We assume that UHECRs are protons injected with a power-law spectrum by 
extragalactic sources. The injection spectrum in taken to be of the form: 

F{E) cx B-t exp(-£/£ max ), (1) 
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where 7 is the spectral index and E max is the maximum injection energy at the source. 
We fix 7 = 2.6 and -E max = 10 215 eV since these values reproduce well the experimental 
results in the lower energy high statistics region around ~ 10 18 5 eV up to ~ 10 20 
eV (see, e.g. |21|). We focus on events above 4 x 10 19 eV, which are generated at 
z <C 1, therefore, source evolution is only marginally relevant 12 . We simulate the 
propagation of protons from the source to the observer by including the photo-pion 
production, pair production, and adiabatic energy losses due to the expansion of the 
Universe. We assume the Universe has matter and dark energy densities as fractions 
of the critical density given by fl m = 0.3 and Oa = 0.7, respectively. 

Our source distribution in space is generated by a random placement of 
sources for a given spatial density and particles are emitted from randomly chosen 
source positions. The source redshifts are generated with a probability distribution 
proportional to 

dn . . , dt 

dl°^dl< (2) 
where dt/dz gives the relation between time and redshift (see, e.g., the expression 
given in and r(z) is defined as 

r = c / ^77T > ( 3 ) 



R(t) ' 

where t s is the age of the Universe when the event was generated, to is the present 
age of the Universe, and R(t) is the scale factor of the Universe. Once a particle has 
been generated, the code propagates it to the detector calculating energy losses and 
taking into account the detector energy and angular resolution. The source angular 
coordinates on the celestial sphere are chosen randomly from a uniform distribution 
in right ascension and with a declination distribution proportional to cos (54 Since we 
neglect the effects of the magnetic fields on the propagation, we ignore sources outside 
the experimental field of view and assign to visible sources a probability proportional 
to: 

r(z)-M<5), (4) 

where r(z)~ 2 takes into account the distance dependence of the solid angle and oj(6) 
is the relative exposure of the experiment in the given direction. 

For a given simulated set of events, we calculate the angular two point correlation 
function as defined in [191 I23j to study the expected departure from an isotropic 
distributions in the sky. For each event, the number of events within a circle at an 
angular distance 9 and a bin width A9 is summed: 

and divided by the area S(9) — 2ir\ cos(#) — cos(# + A0)| of the angular bin between 
9 and 9 + A9, and Rij(9) — 0, 1 counts the events in the same bin. 

For the Auger observatory in the Southern Hemisphere, we assume a total 
acceptance of 7,000 km 2 sr independent of energy above 10 19 eV. For the angular 
resolution we used 1°, for the energy resolution 20%, and for the exposure dependence 
on the arrival direction we used the analytical estimate given in Ref. 3 . By 

X The declination 6 is denned as being 0° on the equatorial plane and ±90° at the poles, so 
df2 = cos <5d5da. 
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Table 1. Number of events expected with E > 10 eV. 



p (Mpc~ 3 ) 


"20 / lyr 

"19.6 = 70 


"20 / 5yr 
n w , 6 = 350 


"20 / 10yr 
"19.6 = 700 


"20 / 15yr 
"19.6 = 1050 


cont. 


2.9 ± 1.7 


14 ± 4 


29 ± 5 


43 ± 6 


1CT 3 


2.5 ± 1.3 


13 ±3 


26 ±6 


39 ±8 


lcr 4 


2.0 ± 1.4 


10 ±4 


22 ±6 


34 ±7 


icr 5 


1.6 ± 1.5 


7.3 ±3.3 


14 ±5 


20 ±7 


10~ 6 


0.4 ±0.7 


2.1 ± 1.7 


4.2 ±2.9 


7.2 ±4.3 



comparison, AGASA had an acceptance of 160 km 2 sr and reported 886 events above 
10 19 eV over a decade with an accumulated exposure of 1, 645 km 2 sr yr. If the two 
experiments had the same energy calibration, Auger would accumulate over ~ 10 4 
events above 10 19 eV in 5 years of operation. However, the first released Auger 
spectrum [21] suggests that the energy calibration of AGASA is systematically 
high and, consequently, a lower flux is observed for a fixed energy at Auger. 

We normalize our simulations to the Auger flux above 4 x 10 19 eV |23] (indicated 
in Table^as nig. 6: number of events above 4 x 10 19 eV) and determine the number 
of events with energy above 10 20 eV as a result of our simulations. These numbers 
are summarized in Tabled f° r a continuous distribution of sources and for different 
values of the source density. 

The numbers in Table ^ require some further comments: the spectrum of cosmic 
rays as detected by Auger at the present time is rather odd, in that it appears to 
have a dip instead of a bump at energies around 4 x 10 19 eV. This may easily lead 
to an incorrect estimate of the number of events expected at higher energies. In Fig. 
^we show the Auger data and the predicted spectrum as normalized at low energies 
(10 19 eV, dashed line) and at 4 x 10 19 eV (solid line). Systematic uncertainties on the 
flux and on the energy determination are indicated by double arrows at two different 
energies |2"4] . 

The figure clearly illustrates the fact that the normalization at 4 x 10 19 eV may 
cause an underestimate of the number of events expected at higher energies, since 
most likely the feature observed at that energy is the result of a statistical fluctuation. 
In order to keep memory of this problem we carry out our analyses for both instances 
of normalization. 

The number of events as obtained with the higher normalization are reported in 
Table an are systematically higher by a factor ~ 2. 

We parametrize the mean density for the distribution of sources by p — 10 — r 
Mpc~ 3 , and choose r between 3 and 6. The mean densities of extragalactic 
astrophysical accelerators of UHECRs should be well represented by this range which 
covers galaxies, groups and clusters of galaxies. For example, black holes in centers 
of normal galaxies should have r ~ 3 while those in active galaxies should have 
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Figure 1. Data points from the Pierre Auger Observatory 1241 compared with 
theoretical predictions normalized at different energies. 



Table 2. Number of events expected with E > 10 eV using the normalization 
of the spectrum at low energy. 



p (Mpc~ 3 ) 


"20 / lyr 
nig.e = 120 


"20 / 5yr 
n w , 6 = 600 


"20 / 10yr 
nig.e = 1200 


7i20 / 15yr 
7719.6 = 1800 


cont. 


5.0 ±2.2 


25 ±5 


50 ±7 


75 ±8 


1(T 3 


4.2 ± 1.8 


22 ±5 


44 ±9 


67 ± 11 


i(r 4 


3.6 ±2.0 


18 ±5 


38 ±9 


58 ± 12 


10" 5 


2.5 ± 1.8 


13 ±5 


24 ±8 


35 ± 10 


icr 6 


0.8 ± 1.0 


3.5 ±2.3 


7.4 ±4.2 


13 ±7 



r > 6. Similar densities are expected for rich clusters of galaxies. The most common 
extragalactic objects that may house UHECR sources are galaxies, which have p 
depending on their type and luminosity ranging from r = 2 to 3. For example, 
the Sloan Digital Sky survey recently reported a comprehensive study of galaxy 
clustering in a large volume limited sample at relatively low redshifts (z between 
0.015 and 0.1). In their samples, galaxy number densities vary between 2 x 10 h 3 
Mpc~ 3 (in their sample with limiting r-mag nitude -18) to 6 x 10~ 3 h 3 Mpc -3 (for the 
brighter sample). The number density in groups of galaxies range from 6 x 10 h 3 
Mpc~ 3 (starting with groups of 3 galaxies) to richer clusters with 10~ 7 h 3 Mpc -3 . 
Choosing the Hubble parameter, h = i?o/100 km/s/Mpc = 0.75, the observed p range 
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Figure 2. Two point correlation function expected for Auger events above 
4x 10 19 eV after 1, 5, 10, and 15 years, for continuous distribution of sources (black 
circles), and number densities of 10 — 3 Mpc -3 (cyan stars), 10 — 4 Mpc — 3 (green 
downward triangles), 10 -5 Mpc -3 (blue squares), and 10 — 6 Mpc -3 (red upward 
triangles). No limitation on the minimum distance for sources were imposed. 
Error bars are asymmetric 1 a. The data points in this and in all subsequent 
figures have been slightly offset (horizontally) for clarity of presentation. 



from 8 x 10~ 3 Mpc" 3 to 4 x 10~ 8 Mpc -3 . As we discuss below, detecting SSAs for 
p > 10~ 3 Mpc~ 3 will require many years of observations, while p <C 10~ 6 Mpc -3 may 
generate higher UHECR clustering than reported so far, unless magnetic fields move 
the threshold for charged particle astronomy to energies around 10 20 eV. Such low 
densities of sources however would generate a sharp GZK cutoff in the diffuse spectrum, 
starting at energies even below 10 20 eV. This situation appears to be disfavored even 
on the basis of current data. 



2.1. Two point correlation function above 4 x 10 19 eV 

Figure |21 shows the two point correlation function of simulated events with energies 
above 4 x 10 19 eV for the full Auger South aperture and exposures after 1, 5, 10, and 
15 years of operation. The number of events is normalized at the Auger data at energy 
4 x 10 19 eV. The black circles show a continuous distribution of sources, while sources 
with number densities of 10 -3 Mpc -3 are represented by cyan stars, 10 -4 Mpc -3 
by green downward triangles, 1CT 5 Mpc~ 3 by blue squares, and 10~ 6 Mpc -3 by red 
upward triangles. The error bars are asymmetric 1 a standard deviations calculated in 
each energy bin for events above and below the mean. It is clear from this figure that 
after one year with the full Auger South aperture a discrete distribution with p < 10~ 4 
Mpc -3 should be distinguishable from a continuous source distribution. After 5 years, 
even ap~ 10 -3 Mpc -3 can be detected at the few a level. After a decade, departure 
from continuous should be detected at many tens of a. 

Although departures from continuous distributions should become clear soon, 
distinguishing between different specific source densities in Fig. [21 is more difficult. 

In the 100 realizations for the assumptions in Fig. there are often configurations 
with sources within 10 Mpc. These randomly placed nearby sources produce a 
significant dispersion in the two point correlation function at small opening angles 
and generate an overlap between different choices of p. To help distinguish between 
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Figure 3. Same as Figl^lbut the minimum distance for nearby sources is rf m i n (p)- 



different number densities, we impose a lower limit in the distance to the nearest 
source and show the constrained angular correlations in Fig. In Fig. |2| the nearest 
source is at a distance larger than d m i n (p), where 

dmin(p) = ~ p- 1/3 • (6) 

As can be seen in Fig. |2l imposing a minimum distance alleviates the degeneracies 
between different source densities after 5 years of full Auger South operations. After 
the first full year exposure, number densities above 10~ 4 Mpc~ 3 give clear departure 
from continuous distributions, while for 10 -3 Mpc -3 , 5 years are necessary. 

In Table|2| we list the mean number of doublets (Nd) expected after 1, 5, 10, and 
15 years of Auger South operations for E > 4 x 10 19 eV and for r between 3 and 6. In 
the same table we report, along with TVd, the asymmetric la errorbars for realizations 
with a number of doublets above or below the mean. The opening angle used to count 
doublets is 1° and the realizations have minimum source distance d m i n (p) as in Fig. 
121 It is worth stressing that jVa in the table is the total number of doublets, including 
doublets inside possible higher multiplicity clusters. 

If the number of events used in the simulation is normalized to the low energy 
spectrum as measured by Auger South, then the number of doublets expected in the 
future Auger operation are those reported in Tabled 



2.2. Two point correlation function above \Q 20 eV 

UHECR clustering and source positions should become easier to identify with events of 
energies > 10 20 eV, given that at these energies UHECRs are most likely protons and 
their trajectories are less likely to be affected by extragalactic and galactic magnetic 
fields than lower energy events. The difficulty is clearly in the limited statistics of 
events that can be accumulated in this energy range by Auger South, since the flux 
is a steeply decreasing power of energy and energy losses that give rise to the GZK 
feature are also at play in this energy range. A larger array as currently being discussed 
for the Auger North Observatory should help significantly the ability to do charged 
particle astronomy at this energy scale. 

We simulated the expected angular correlation function in Auger South for events 
with energies above 10 20 eV for different source densities. In Fig. the correlation 
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Table 3. Number of doublets for E > 4 X 10 19 eV, using the normalization of 
the number of simulated events to the number of events observed by Auger above 
4 X 10 19 eV. 



p (Mpcr 3 ) 


^Vd / lyr 


N d 1 5yr 


iV d / 10yr 


7V d / 15yr 


cont. 


33 +0 ' 77 






P- ,+8-1 
67 ' 5 -8.3 


1(T 3 


1.07 +2 ' 75 

-0.80 








1(T 4 




™™ 






10" 5 








2 ' 181 -530 


1(T 6 








9 946 +4 ' 681 

J ' y -2,010 



Table 4. Number of doublets for E > 4 X 10 19 eV, using the normalization of 
the number of simulated events to the number of events observed by Auger at low 
energy. 



p (Mpc" 3 ) 


^Vd / lyr 


Wj / 5yr 


N d 1 lOyr 


«d / 15yr 


cont. 


92 + °" 

U ' -0.92 






200™ 


1CT 3 


3 - 36 -2.03 




366 -109 




io- 4 






i.»i!S 




1CT 5 










1(T 6 


114- 









function is shown for sources located at minimum distances (i m i n (p) as in Fig. [3] In 
addition, we fixed the number of events above 10 20 eV to the mean number given 
in Table ^ If we let the number of events fluctuate between realizations, the error 
bars are significantly enhanced. In principle, once Auger has run for 5 years, we will 
know the number of events above 10 20 eV. As discussed below, the number of events 
expected above 10 20 eV for a fixed injection spectrum, depends on the source density. 
Combining a study of the GZK feature with the small scale anisotropics should yield 
p. As shown in Fig. |3| a departure from a continuous source distribution is detectable 
at the 1 a level after 5 years for p < 10 -5 Mpc -3 . Larger densities require 10 to 15 
years before clear detectability. Even after a decade, distinguishing between different 
source densities will be a major challenge. 

In Table El we show the predicted number of doublets at energies above 10 20 eV 
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Figure 4. Two point correlation function for events above 10 20 eV after 5, 10, and 
15 years at Auger South for a continuous distribution of sources (black circles), 
and number densities of 10 -3 Mpc -3 (cyan stars), 10 — 4 Mpc -3 (green downward 
triangles), 10~ 5 Mpc -3 (blue squares), and 10 -6 Mpc -3 (red upward triangles). 
The minimum distance for the nearest source is rf m i n (p). 



Table 5. Number of doublets for E > 10 20 eV with normalization to the Auger 
data above 4 X 10 19 eV. 



p (Mpc~ 3 ) 


N d / 5yr 


N d 1 lOyr 


N d 1 15yr 


cont. 




n 0.+ 1 ' 38 
- 04 -0.04 




10~ 3 




3 05 +4 ' 84 




10~ 4 


2 32 +3 ' 2S 
Z * -1.S0 


11.1 + 137 
-6.1 


24.S+ 23 " 

-12.6 


10~ 5 






44.4 

-19.4 


10~ 6 






+94 
125 -27 



after 5, 10, and 15 years of operation of Auger south for different values of the source 
density. 

As stressed in the previous section, the number of events at energies above 10 20 
eV is estimated by using as a normalization the flux currently observed by Auger 
South at 4 x 10 19 eV. However, as shown in Fig. ^ the observed spectrum appears 
to be at odds with the observed low energy spectrum. If the solid line in Fig. \I\is 
used for the normalization of the number of events above 10 20 eV, then the situation 
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Figure 5. Two point correlation function for events above 10 20 eV after 5, 10, and 
15 years at Auger South for a continuous distribution of sources (black circles), 
and number densities of 10 -3 Mpc -3 (cyan stars), 10 — 4 Mpc -3 (green downward 
triangles), 10~ 5 Mpc -3 (blue squares), and 10 -6 Mpc -3 (red upward triangles). 
The minimum distance for the nearest source is e£ m i n (p). 



Table 6. Number of doublets for E > 10 2t) eV with normalization to the low 
energy Auger data. 



p (Mpc" 3 ) 


N d / 5yr 


N d / 10yr 


N d / 15yr 


cont. 


02 + °- 98 

U - -0.02 




39 +1 ° 5 

U '^ y -0.39 


10~ 3 


3 29 +5 ° 8 






1(T 4 


10 - 5+ 53 
— 5.3 






icr 5 


23 9 +33 4 






icr 6 


59.6_ 17g 


228!™ 


503!- 



improves somewhat (the number of events above 10 20 eV is roughly doubled in this 
case). This is shown in Fig. where we plot the two point correlation function for 
the new number of events above 10 20 eV after 5, 10, and 15 years. 

The number of doublets obtained for this second instance of normalization to the 
data is shown in Table 
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Figure 6. Energy spectrum for different source densities compared with Auger 
data (open squares). The sources densities are f0 — 6 Mpc -3 (red upward 
triangles); 10 -5 Mpc -3 (blue squares), 10 -4 Mpc -3 (green downward triangles); 
f0~ 3 Mpc -3 (cyan stars); and a continuous distribution of sources (black circles). 



2. 3. Spectra for different number densities 

As shown in |Ki| the combination of small angle clustering and spectral informations 
is a powerful tool to investigate the nature of the sources of UHECRs. A low density 
of sources would in fact be responsible for more evident small angle clustering of 
the arrival directions, but at the same time would also induce a sharper shape of 
the GZK feature, and therefore a lower number of events at the highest energies. A 
larger density of sources has the opposite effect: the GZK feature is shallower (more 
events at the highest energies) but the clustering is mild. As pointed out in the 
AGASA data on the small scale anisotropics appear to be at odds with the spectrum, 
as measured by the same experiment. 

Here we investigate the potential to use a similar technique for the data of the 
Pierre Auger Observatory. 

In Fig. we show the energy spectrum around the GZK feature that may be 
detected by Auger South after 15 years of full aperture operations, assuming the 
normalization at low energy. The spectrum is shown for different source densities 
ranging from a continuous source distribution (black circles) to source densities 
between 10 -3 Mpc -3 (cyan stars) and 10~ 6 Mpc -3 (red upward triangles). Also 
shown is the first release of the Auger spectrum (open magenta squares) |24j and an 
analytical prediction (continuous line). 

The number of events at energy above 10 20 eV is a rather sensitive function of 
the source density, as can also be understood from Tables and |21 For the case of 
normalization at low energy, one could have as many as ~ 75 events above 10 20 eV 
in 15 years of operation of Auger South, for a continuous (and unlikely) distribution 



Small Scale Anisotropy Predictions for the Auger Observatory 



13 



of sources, or as little as ~ 13 events in the same energy region if the sources have 
a mean local density 10 -6 Mpc -3 . Unfortunately this also means that the statistical 
uncertainties in the fluxes at the highest energies become large if the source density 
is too low, and it becomes correspondingly harder to have a precision measurement of 
the shape of the GZK feature in these cases. For source densities between 10 _6 Mpc~ 3 
and 10 _5 Mpc- 3 Au ger is expected to become statistics dominated at energies ~ 10 
eV. However at this energy the different predicted spectra already differ from that 
generated by a continuous source distribution by about a factor ~ 10. Again, the 
combination of spectra measurements and anisotropy measurements can play an 
instrumental role in inferring hints to the nature of the sources of UHECRs. 

3. Conclusions 

The imminent completion of the Auger Observatory in Argentina will mark the 
beginning of a new era in cosmic ray astrophysics. The combination of a very large 
ground array and fluorescence telescopes will provide a detailed measurement of the 
spectrum at the highest energies. 

One of the main open questions in Cosmic Ray Physics is the presence of a GZK 
feature in the spectrum. If this feature is in fact there, it becomes very important to 
measure its shape, which can provide us with information on the number of sources 
and their spatial distribution. Independently of the presence or absence of the GZK 
feature, the problem of finding an acceleration mechanism and a class of sources that 
may harbor the accelerator remains of paramount importance. Auger will most likely 
clarify whether the spectrum of UHECRs has or not the GZK suppression. In order 
to answer the second big question, we need high quality information (namely high 
statistics) about the directions of arrival of UHECRs, and their chemical composition. 
The first hints of such new era of charged particle astronomy are likely to come from 
small deviations from isotropic distributions in the next few years of Auger South 
operations. 

In this paper we made quantitative predictions of the amount of small scale 
anisotropies expected in Auger in the next decade for a range of plausible source 
densities. We investigated two instances of normalization of the number of events 
expected at the highest energies, depending on whether the procedure is applied 
to the data measured by Auger at 4 x 10 19 eV or rather at lower energies where 
statistical errors should be less important. The amount of small angle clustering has 
been quantified through the two point correlation function, both for events above 
4 x 10 19 eV and above 10 20 eV. We conclude that it will be sufficient to run Auger 
for a few years to detect the departure of the two point correlation function from that 
expected from a continuous distribution of sources. On the other hand, it will be more 
problematic to measure the mean source density with high accuracy, because of the 
relatively low number of events and correspondingly large statistical fluctuations. As 
already proposed in [T3], the combination of small angle clustering and spectrum of 
diffuse cosmic rays is able to provide us with more stringent constraints on the source 
density. This will be true for Auger as well, although the statistical error bars on the 
spectrum for the case of discrete sources, as shown in Fig. are expected to become 
relatively large at energies around 10 20 ' 2 eV. This result is rather more pessimistic 
than the ones previously obtained in because such previous results were obtained 
adopting the normalization (not the shape) of the AGASA spectrum, which has now 
been found to lead to cosmic ray fluxes about three times larger than those actually 
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measured by Auger at the present time. 

Despite these problems, it is fair to conclude that once the departures from 
isotropy will be detected, the estimate of the mean density of sources, p will be within 
reach. The even more challenging possibility of measuring the spectrum of a single 
source ^2] w iU require a significant increase in statistics at the highest energies as 
currently being discussed in the context of the Auger North Observatory or future 
generation cosmic ray observatories, that will most likely be located in space. 
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